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LOGARITHMIC  DERIVATIVES  OF  TWO  NORMALIZING  FUNCTIONS 

:  *  Geoffrey  S.  Watson 

Princeton  University 

ABSTRACT 

If  the  unit  vector  x  In  JRq  has  a  probability  density  proportional 
to  exp  k  u'x  ,  or  exp  ic(u»  .  then  the  statistical  theory  for  x  depends 
largely  upon  the  logarithmic  derivatives  Aq(ic)  and  of 

/  eKtU-t2)vdt  and  /  eKt2(l-tZ)vdt 
-1  -1 

where  v«(q-3)/2  .  This  paper  gives  a  self-contained  study  of  the  functions 
Aq(ic)  ,  Bq(ic)  ,  of  the  computational  problems  of  calculating  these  functions,  of 
solving  y-AqU)  and  y-B^tic)  ,  and  of  finding  the  variance  stabilizing  trans¬ 
formations  /'A'UJ^k  ,  are  also  discussed. 


Key  words:  Flsher-von  Mlses  distribution,  Scheldegger-Watson  distribution, 
power  series,  asymptotic  expansions,  variance  stabilizing 
transformations,  Bessel  functions,  Runner  functions,  Rlccatl  equation. 
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1.  INTRODUCTION 


Two  probability  densities  on  the  surface  of  the  unit  ball  In  ®q 


are  given  by 

fj(x)  ■  ajU)"1  exp  k  v'x  (1.1) 

f2(x)  -  bjUf1  exp  k(ii-x)2  (1.2) 

when  y*x  Is  the  scalar  product  of  the  two  unit  vectors  x  and  y  .  y  Is 
called  the  modal  direction  and  k  a  concentration  parameter.  In  (1.1)  k>0  . 

In  (1.2).  k  may  be  any  real  number.  If  ic  +  0  either  density  becomes  uniform. 

The  statistical  theory  of  (1.1).  the  Flsher-von  Mlses  distribution,  has 
been  studied  most  recently  by  Watson  (1981a,  1981b).  The  statistical  theory  of 
(1.2),  the  Scheldegger-Watson  distribution,  has  been  studied  for  general  q  In 
Watson  (1981c).  It  Is  found  that  the  theory  turns  upon 


Aq(ic)  -  a^(K)/aJ(<)  , 


Bq(*0  •  ^(Kj/b*^)  , 


(1.3) 


(1.4) 


respectively.  The  Integral  definitions  of  a*  and  b*  can  be  simplified  by 
writing  x  •  ty  ♦  (l-t2)1^  where  t  -  y'x  ,  and  (  Isa  unit  vector  orthogonal 


to  y  .  Then,  writing 


we  have 


«q  ■  area  of  nq  ■  2irq^2/r(q/2)  , 


•;(«)  *  VlV**  *  b$*K)  "  Vl  bqW 


(1.5) 


A 
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where 

ILlic)  •  /  «'ct(l-t2)(q“3,/2dt  ,  11.6)  * 

7  -1 

b  (k)  -  /  eKt^l 1-t2) tq"3)/2dt  .  (1.7) 

q  -1 

and 

Aq(ic)  -  a'(<)/aq(.c)  ,  Bq(<)  »  b'(x)/bqU)  (1.8) 

The  distribution  (1.1)  originated  in  Statistical  Physics  with  Langevin's 
work  on  magnetism.  (Our  nomenclature  refers  only  to  the  first  uses  of  (1.1)  for 
statistical  Inference.)  Dyson,  Lleb  and  Simon  (1978)  In  proving  the  existence 
of  spontaneous  magnetization  at  sufficiently  low  temperatures  needed  some  properties 
of  log  aq(ic)  .  Suppose  h(y)  Is  defined  by 

h(y)  *  log  /  exp  y'x  dy(x) 

“q 

where  y(x)  Is  any  measure  on  nq  .  Then  h(y)  Is  a  convex  function  of  y  . 

If  c  Is  the  largest  eigen  value  of  C32h/3y,j3y  3  ,  -y  cy'y  ♦  h(y)  is  concave 
In  y  .  If  y  Is  the  uniform  measure  on  nq  .  c-q"1  ,  a  proof  of  which  Is  given 
below.  Their  first  result  Is  basic  for  Laplace  transforms  —  see  e.g.  Bamdorff- 
Nlelsen  (1978). 

Schou  (1978)  gave  seme  properties  of  and  expansions  for  Aq(ic)  by 
recognizing  that  aq(x)  is  proportional  *(q/2)-i  »  (k)k’^2^’1  and  using  the 
properties  of  the  modified  Bessel  function  of  the  first  kind,  Iy(x)  .  Watson 
(1956)  gave  some  results  on  Bq(ie)  ,  q»3  . 


He  wlU  derive  all  the  properties  of  Aq(ic)  (Section  2)  and  Bq(ic) 
(Section  3)  directly,  rather  than  as  the  ratios  a'/a  ,  b'/b  .  More  properties 
and  more  terns  In  the  expansions  will  be  given  than  heretofore.  Furthermore 


In  practice  one  needs  to  be  able  to  compute,  any  q  and  k  ,  Aq(<)  ,  Bq(ic)  , 

solve  the  equations  y  «  Aq(ic)  ,  y  ■  Bq(<)  and  to  apply  the  variance  stabilizing 

transformations  (derived  In  Watson  1981a,  1981b)  gq(tc)  ■  /  A'(ic)^dk  . 

hqU)  »  /  B'(ic)^dk  numerically.  These  matters  are  discussed  In  Section  4. 

While  we  have  wished  to  emphasize  that  It  Is  Aq  and  Bq  rather  than 

8^  and  bq  which  matter,  we  have  for  completeness  Included  all  known  properties 

of  the  functions  aM  and  b.  .  The  former  Is  associated  with  the  modified 

<1  Q 

Bessel  function  of  the  first  kind  Iv(z)  and  the  latter  to  the  less  well  known 
Kumner  function  of  the  first  kind  M(a,b,z)  . 
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2.  The  functions  aq(ic>  .  Aq(ic)  and  Its  Inverse. 
Consider  the  functions  of  k>0  . 


a.(ic)  -  /  eKt(l-t2)tq*3,/2dt  ,  q>2  , 
q  -1 

(2.1) 

AqU)  *  a'(ic)/aq(<)  . 

(2.2) 

For  the  applications  now  envisaged  q  Is  Integer  >2 

but  the  following  results 

are  true  for  real  q>2  . 

It  Is  clear  that  If  q^ > q2  *  0<a  U)<a  (k) 

.  It  Is  simpler  to  work 

with  1 

yv(tc)  ■  aq(ic)  ,  v  •  (q-3)/2>-%  , 

(2.3) 

Mv(k)  •  y'(ic)/yv(K)  ■  Aq(ic) 

(2.4) 

Examination  of 

y'{<)  -  /  eKt(l-t2)vdt 
v  -1 

shows  that 

o<y'(x)<yv(K) 

(2 .3) 

whence 

0<Aq(w)  -  Hv(w)<l 

(2.6) 

Elementary  manipulations  show  that 


K  '  7T&TT  Vl  •  <v>-»  • 

! 

f 

I 


(2.7) 


«W»  . . 


\ 
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.  syj  ■  -(2v+l)yv  ♦  2v  yv-1  ,  (v>0)  * 

i- 


(2.8) 


•  ny''*  2(v+l)y'  -  icyv  -  0  ,  (v>-%)  .  (2.9) 

It  «ay  be  verified  that  (2.9)  has  a  solution  proportional  to  k“^v+>^ 

whence 

aq(ic)  «  (2ir)q/2  Iq/2-i<*c>  *’q/2+1  (2.10) 


Here  Iv(ic)  Is  the  modified  Bessel  function  of  the  first  kind  as  defined  In 
Matson  (1952).  Thus  as  k-*0  , 

«  (k>  *  (2n)q/2  k^2*1  £  (K/2)q/2‘1+2r  (2.11) 

q  r-0  ntq)nvfTr 

while  as  , 

aq(ic)  •  (2vr)q/2'1  K(q_3,/2  eK(l  +  Ofic"1))  (2.12) 

Since 

A'<*)  .  -  («'/«q)2  .  (2.13) 

S'-  *q  • 

It  follows  that 

A'(k)  ♦  A2(k)  <  1  .  (2.14) 

Using  the  Cauchy  Inequality*  ajj'aq>(aq)2  ,  In  (2.13)  and  the  result  of  Qyson* 
Lleb  and  Simon  (1978)*  we  have 


0<A£(.c)<l/q 


(2.15) 
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Returntng  the  probabilistic  background,  t-y'x  where  x  has  the 
distribution  (1)  so  that 


Aq(0  »  Et  . 

(2.16) 

A'Oc)  •  Et2  -  (Et)2  »  var  t  , 

(2.17) 

so  that(2.6)  and  (2.15)  are  statistically  obvious.  For  Aq(ic)  »  var  t  will  be  a 

maximum  when  k*0  because  of  (2.18)  below.  It  Is  then  equal  to  1/q  by  a  direct 

easy  calculation.  Further  since 

a"a'  a'  3 
»  -J—  -  3  -V*  +  2(-S.)  . 

q  q  aq  q 

*  Et3  -  3  Et2Et  +  2(Et)3  , 

-  E(t-Et)3  . 

so  from  the  evident  skewness  of  the  distribution  of  t  for  k >  0  it  follows  that 

A~(k)  <  0  .  (2.18) 

Thus  the  function  Aq(ic)  1$  non-decreasing  and  convex  on  (0,®)  taking  Its 
minimum  at  k«0  and  maximum  as  ie«®  In  the  range  CO, 13  .  While  we  will  not 
went  Aq(ie)  for  negative  k  ,  it  is  helpful  to  observe  that  Aq(ic)  Is  an  odd 
function  of  k  .  It  Is  easy  to  show  that  Aq.A'.A''  ....  are  the  successive 
cumulants  of  t"u'x  ,  a  fact  which  Is  used  in  Watson  (1981d)  to  derive  Edgeworth 
expansions  for  certain  limiting  distributions. 

Integrating  a'*  by  parts  and  using  (2.3).  we  find  that  Aq(ic)  satisfies 
the  Rlccatl  equation 

9 
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A'(k)  ■  1  -  Aq(ie)  -  Sjl  Aq(ic)  .  (2.19) 

as  found  via  Bessel  functions  by  Schou  (1978).  Similar  manipulations  show  that 

v*>  •  ¥  •  f"*3)  •  ‘2-3<» 

as  In  Schou  (1978) (Al) .  By  (2.6),  (2.15),  and  (2.8)  It  follows  that 

A'(-)  *  0  ,  Aq(»)  ■  1  (2.21) 

Inspection  of  AqU)  as  k*0  shows  Aq(0)  ■  0  so 

Aq(»)  »  £  .  Aq(0)  -  0  ( 2.22 ) 


He  now  give  more  detailed  expansions  for  AqU)  as  k  +  O  and  by  solving 

(2.19).  The  difference  equation  will  be  used  later  to  suggest  how  to  tabulate 
Aq(ic)  for  various  values  of  q  . 

For  q«3  ,  there  Is  a  simple  explicit  formula 

Aj(k)  ■  coth  tc  -  k"1  (2.23) 

For  q«2  ,  we  may  write 

AgU)  -  IjUVIqOc)  (2.24) 

since  generally 

Aq(«e)  -  Iq/Z+l^^q^-lto  (2-24') 

For  q*2  ,  the  zeros  of  Iq  nearest  to  the  origin  are  at  ±2.411  (approx.).  For 

q*3  ,  the  zeros  of  I^2  nearest  to  the  origin  are  at  ±3.811  (approx.).  These 


•8- 


comfits  help  to  define  the  radius  of  convergence  of  the  series  expansion  (2.25) 

below. 

For  it*0  ,  we  way  substitute  ajie  ♦  a3ie3  ♦  agic5  +...  for  the  (odd) 
function  Aq(ic)  In  (2.19).  We  then  find  that 


A_(ic)  *  “  <  -  — W— K3  + 
q  q  qZ(q+2) 


q4(q+2)(q+4) 


k5  ♦  0(k7) 


(2.25) 


where  the  first  two  terms  agree  with  Schou  (1978). 

For  K+~  ,  we  may  substitute  a  series  In  powers  of  ic"1  In  (2.19)  and 
find  that 

+  4  +  0(4)  (2.26) 

where  the  first  three  terms  agree  with  Schou  (1978).  For  q«3  the  terms  In  *"2 
-3 

and  k  are  zero.  But  from  (2.23),  we  may  write 


Aj(ic)  "  (1  -  ^  +  “  e"2,c)(l  -  e"2lc)"1  (2.27) 

so  we  see  that  A3(k)  -  1  +  £  Is  exponentially  small. 

From  the  form  of  (1)  It  Is  obvious  that  maximum  likelihood  leads  to 
equations  for  k  of  the  form  y  »  A(k)  so  that  we  need  k  •  A“*(y)  for  all 
y  In  C0,13  .  The  equation  (2.19)  can  be  rewritten 

<(l-y*)K  -  (q-I)y)  fy  *  *  (2.28) 

3 

Since  the  Inverse  function  will  also  be  odd  we  may  put  k  •  a^y  +  a^y  ♦... 

In  (2.28)  for  y  small.  He  find  then,  as  y*0  , 
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*  ■  «>  y +  y3  *  rq+2^1) y5  +  0(y7) 


(2.29) 


Alternatively  Me  could  revert  the  series  (2.25)  to  find  (2.29).  To  find  ic  when 
y  Is  near  unity  we  may  revert  (2.26)  using  z  »  1-y  «  1-Aq(ic)  or  set  z  »  1-y 
In  (2.19)  and  solve  for  a  series  In  Inverse  powers  of  z  .  Either  way  we  find 
that  the  solution  of  y  «  Aq(ic)  for  y  near  unity  Is  defined  by 


-  -  qlT  (*-y)  +  Tq^fy  (1-y)2  +  {1"y)3  +  0{(1‘y)4)  (2*30) 


The  case  q*3  Is  again  special  but  (2.27)  yields 


J  ^(l-y)  +  2e"'c  , 


J  ~ (1-y)  +  2  exp  -  (1-y)"1  . 


(2.31) 
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3.  The  functions  bq(K)  ,  Bq(K)  and  Its  Inverse. 
Here  we  consider  the  functions 

t 

bqU)  -  zv(k)  •  /  eKt  (l-t2)vdt  , 


Bq<*>  «  *  bj/bq  -  z-/zv  . 


13.1) 


13.2) 


where  v»(q-3)/2>-%  and  —  <k<»  . 

The  functions  zv(»c)  ,  bq(<)  and  their  successive  derivatives  form 
descending  sequences  because  |t|<l  .  rurther 


0<b<I1  <  %  1f  1 


0<z  <  z  If  v,  > 


V1  v2 


l>q2 


0  <  Nv(ic)  -  Bq(<)  <  1 


Elementary  calculations  show  that 


13.3) 


*  *v  -  Vi  •  (v>-%)  . 

2k  z '  -  2v  zv-l  -  (2v+l)zv  ,  (v>0) 


2k  z"  *  (2v+3  -  2k)z'  -  zv  -  0  ,  (v>-%)  . 

In  this  case  we  will  proceed  without  identifying  and  further  describing  the 
function  zv(k)  ■  bqU)  .  Details  are  given  in  the  Appendix. 


(3.4) 


(3.5) 


(3.6) 
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From  Cauchy's  Inequality  and 

B'U)  •  b"/bq  -  (b'^J2  .  (3.7) 

b''(ie)  <  bq(<)  , 

It  folloMS  that 

0  <  B'U)  <  1  -  B2(k)  13.8) 

2 

If  the  random  vector  x  has  the  distribution  (2)  and  u  *  (ji'x)  » 


then  1  >  Eu  >  Eu2  >. . .  .>  0  and 

Bq0c)  -  Eu  .  (3.9) 

BqU)  ■  Eu2  -  (Eu)2  *  var  u  (3.10) 

B~U)  -  Eu3  -  Eu2Eu  -  2Eu(Eu2-(Eu)2)  -  E(u-Eu)3  .  (3.11) 

From  the  skewness  of  the  distribution  of  u  It  follows  that 

B~(ic)  >  0  (ie<0)  Bq'(ic)  <  0  (k>0)  (3.12) 


Thus  Bq(ic)  Is  non-decreasing  and  concave  on  (-*,0)  ,  convex  on  (0,»)  ,  taking 
Its  minimum  at  and  Its  maximum  at  +»  In  the  range  CO, ID  . 

If  we  put  tc*0  In  (3.4),  and  In  (3.5)  with  v+l  Instead  of  v  ,  we  find 

that 


Bq(0)  -  1/q  . 


(3.13) 
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Then  (3.4)  shows  that 

$q(0)  »  2/q(q+2)  .  (3.14) 

From  (3.4)  and  (3.S)  we  find  the  recurrence  relation 

Bq(«)U-Bq.2(K))  -  ♦  ^fBq-2(ic)  (3.15) 

> 

Now  Bq(ic)  and  Bq-2(ic)  tend  to  limits  (the  same  from  (3.4))  In  ,  13  as 
ic which  (3.15)  shows  to  be 

Bq(")  “  1  (3.16) 

Similarly  Bq(ic)  lies  In  CO  ,  ^3  for  *<0  and 

Bq(-)  -  0  (3.17) 

Finally  we  note  that  Bq(ic)  satisfies  the  Rlccatl  equation 

1  -  2*c  B-(<)  -  2k  B*(ic)  -  (q-2x)Bq(,c)  (3.18) 

which  holds  for  q>2  . 

To  find  expansions  for  Bq(<)  as  io*-0  ,  ,  and  ,  we  may 

here  use  either  the  difference  equation  (3.15)  or  the  differential  equation  (3.18). 
He  know  of  no  explicit  forms  for  special  values  of  q  but  we  will  discuss  the 
cases  q*2  and  q»3  separately.  After  harder  calculations  than  the  last 
section  we  find: 

for  ic-*0 

B-(k)  .  1 .  -SiaiU  K  .  .  (  zSIa-ul. 

q  q  qZ(q*2)  q3(q+2)(q*4)  q4(q+2)(q+6) 

0 

+  8(q-;)Lq"?l2  J  k3  4  0(/)  f  (3.19) 

q4(q+2)(q*4)(q+6) 


13- 


f  or  k -*■+«• 


BqU)  •  1  -  i  -  V-  7  -  js+H-fy  • 


7 


(3.20) 


for  k  ♦  -• 


(K).4i.la^l  .  in 3^41  4+0(-V) 

IC  1C  1C 


Vc,“"7 


(3.21) 


Finally  we  need  expansions  for  the  solution  of  x*Bq(y)  for  y  near 
zero,  near  1/q  and  near  unity.  We  may  revert  the  appropriate  expansions  above 
or  rewrite  (3.18)  as 

{1  -  2*y2  +  (2k:  -  q)y>  ^  »  2k  (3.22) 


For  a  solution  when  y  Is  near  q"1  ,  the  needs  to  be  changed  to 

z  ■  y-q"1  .  We  find 


(3.23) 


r 

t 


As  y-*-l  ,  ic-*-®  and  we  have 


I 

K 


2 

(q-l)Z 


d-y)2 


-  ^4  (1-y)3  +  0((l-y)4) 

(q-1)3 

As  y-*-0  ,  k^-*  and  we  find 

.  I  .  2y  ♦  2(q-3)y2  -  L<l-3)(Sqd9J  y3  «,  0(y4) 

The  formulas  (3.21),  (3.26)  are  unhelpful  when  q*3  . 


(3.24) 


(3.25) 


I 
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For  q«2  ,  the  density  (2)  Is  proportional  to  exp  k  cos29  and  has 
2  2 

normalizing  constant  /  exp  k  cos  9  de  when  k>0  «  the  modes  are  at  9-0,ir 

0 

while  If  ir<b  «  they  are  at  ir/2  ,  3ir/2  .  Thus  there  Is  no  need  to  use 
ic < 0  —  a  rotation  of  ir/2  will  serve.  Further  since  cos20-  |(cos20+l) 

4-26  has  the  density  (1)  with  a  concentration  of  k/2  .  Thus 

B2(k)  -  ^{1  ♦  Ag(£))  (3.26) 

It  Is  easily  verified  that  this  relation  Is  consistent  with  Ag  satisfying 
(2.19)  and  B2  satisfying  (3.18) •  etc..  Thus  we  may  check  (3.19)  against 
(2.25)  and  (3.20)  against  (2.26).  In  so  doing  we  find  that 

B2(k)  "  \  +  f  "  k3  4  0(k5)  (3.27) 

Our  ic3  term  In  (3.19)  was  so  complicated  that  we  hesitated  to  give  It,  but 
this  device  provides  a  check  of  It. 


For  q-3  , 


1  2  ft2  1  rt2 

B,(k)  -  /  re*x  dt  /  /  e*x  dt  . 
3  -1  -1 


If  set  ic--X  ,  X-*-»  .  Then  setting  u  ■  (2X)**t 

2X  B(-X)%1  -  {  l - L_  -i-e’X)1  (3.28) 

fin  (2X)*  fix 

by  noting  the  relation  to  the  Gaussian  and  the  fact  that 

1  -  /*  JL.  e*x2/2  x^„ 

-•  fin  fin 

Thus  (3.28)  explains  (3.21)  l.e.  B3(ie)  ♦  (2k)’1  Is  exponentially  small  as  . 

Similarly  as  ,  w-1  ♦  2y  (see  (3.25))  Is  exponentially  small. 
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♦.  The  computation  of  Aq(ic)  .  Bq(ic)  and  associated  functions. 

From  (2.10)  It  follows  that 

a 

*  *  • 

V*)  "  V2^  /  ^-l^  •  (4*1) 

Since  these  modified  Bessel  functions  are  tabulated  (see  e.g.  Abramowltz  and 
Stegun  (1964)),  A  (k)  may  be  found  for  any  q  and  k  .  To  solve  y»A (k)  we 

”  e  ” 

may  use  Newton's  method  since  (2.19)  , 

A'U)  -  1  -  A*(k)  -  Sji  A^k)  .  (4.2) 

and  (4.1)  enable  us  to  compute  the  derivative.  Thus  hand  calculations  are  not 
a  problem  for  Aq(<)  . 

To  produce  tables  of  the  Aq(ie)  ,  machine  calculations  are  needed.  If 
q  Is  odd,  we  may  first  compute 

Aj(ic)  »  coth  k  -  (4.3) 

and  then  use  the  recurrence  relation  (2.20)  so  that 
Ag(ic)  *  (Aj(ic))  *  -  3tc’*  ,  etc.. 

If  q  Is  odd,  we  need  separate  algorithms  for  A^(k)  and  A4(k)  before  (2.20) 
may  be  used. 

To  compute  A2U)  and  A4(k)  one  may  return  to  the  basic  definition 
and  use  numerical  Integration  to  evaluate  a^ic)  ,  a'U)  for  smaller  values  of 
k  until  the  asymptotic  expansion  (2.26)  comes  In  to  force  and  the  ratio  becomes 
awkward  to  handle.  The  latter  will  not  happen  If  a ^(ie)  and  a'(ic)  are  divided 
by  the  first  terms  of  their  asymptotic  expansions.  From  (2.10)  and  the  expansion 
(Watson  (1952))  of  Iv(k)  •  *«*■•  •  we  have 
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•qU)  (2/«r)<,/2-1  e'  (4.4) 

-  ft 

The  derivative  of  (4.4)  will  be  used  for  a'(ie)  . 

If  Aq(ic)  Is  to  be  tabulated  for  a  specific  q  It  nay  be  best  to 
Integrate  the  differential  equation  (4.2)  numerically.  To  start  the  solution  at 
k«0  •  set 

A^k)  »  tc/q  +  ica(ic)  (4.5) 


so  that  a(ic)  satisfies 


(4.6) 


which  Is  well  behaved  at  the  origin.  Of  course  a(0)»0  .  When  *  becomes  large, 
(4.2)  should  be  used  directly. 

In  Watson  (1981a)  It  Is  shown  that  the  variance  stabilizing  transformation 
for  ^estimates  Is 


9(k)  »  /( A'(k))%dk  .. 


(4.7) 


In  which  the  lower  terminal  Is  arbitrary.  For  q«3  we  have  found  It  convenient 
to  start  at  unity.  The  table  of  g(ic)  Is  then  found  by  numerical  Integration 
using  (4.2)  and  (4.3);  It  1$  more  accurate  for  values  of  k>1  than  for  k<1  . 
Since  k's  <  1  are  rare  for  q>3  In  the  applications  we  have  met,  this  Is  very 
satisfactory.  If  one  regularly  dealt  with  quite  large  k  *  as  In  palaeomagnetlsm. 
It  would  be  satisfactory  to  tabulate  (4.7)  by  numerical  Integration  with  A'(k) 
replaced  by  the  derivative  of  the  asymptotic  expansion  (2.25)  and  starting  at 
k  ■  10  or  20  . 


that 


Turning,  now  to  the  computation  of  Bq(<)  ,  It  Is  shown  In  the  Appendix 


(4.8) 


where  M(a(b,z)  Is  one  of  the  Kunmer  functions.  The  tabulation  most  useful  for 
our  purposes  Is  that  In  Rushton  and  Lang  (1954).  It  gives  both  the  numerator 
and  denominator  of  (4.8)  for  q-1  (1)  5  ,  k  ■  0.02  (•02)*1(*1)1(1)10(10)50,100,200. 
These  tables  may  be  extended  to  negative  ic  by  using  Rummer's  transformation  (A7). 

The  only  known  case  where  (4.8)  reduces  to  more  familiar  functions  Is 
q*3  ,  k<0  .  By  specializing  (4.8)  or  a  direct  calculation. 


u  -4 
e  u  ’ 


du 


(«c>0) 


(4.9) 


and,  setting  for  *<0  ,  -tc*A  , 


b3(x) 


/ 

0 


•u  -4 

■I  * 


du 


U<0) 


(4.10) 


Thus  when  k<0  ,  b3(<)  and  b3(ic)  can  be  found  using  tables  of  the  Incomplete 
ganma- function. 

To  compute  Bq(ic)  for  specific  q^3  and  k  ,  or  to  make  tables.  It 
seems  that  the  best  method  Is  always  to  Integrate  the  differential  equation  (3.18). 
It  may  be  rewritten  as 


BqU)  "  *  £  "  Bq<«>  ’  *  l)  Bq(,c)  l4*U) 

To  start  the  solution  It  Is  essential  to  use  (3.19)  since  (4.11)  Is  badly  behaved 
for  small  k  . 


To  solve  y  ■  BqU)  ,  when  none  of  the  expansions  (3.19),  (3.20)  and 
(3.21)  nay  be  used,  one  must  compute  Bq{  )  and  B'(ie)  from  (4.11). 

To  find  the  variance  stabilizing  transformation 

hU)  •  / (B^(k))1*  dk 

the  range  of  Integration  should  avoid  k«0  If  possible. 


J 
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APPENDIX 


It  Mas  shown  In  Section  3  that  bq(ie)  satisfied  the  differential  equation 
(3.6),  or 


ic  b~U)  ♦  (f  -  k)  b'(<)  -  \  b  (k)  -  0 


(A.1J 


This  Is  Kumner's  differential  equation  —  see  e.g.  Chapter  13  by  L.J.  Slater  In 
Abramowltz  and  Stegun  (1964).  It  Is  self-adjoint  and  a  confluent  form  of  the 
hypergeometric  equation  with  a  regular  singularity  at  ic»0  and  an  Irregular 
singularity  at  .  All  regular  solutions  at  the  origin  of 


zui"  +  (b-z)w'  -  aw  *  0 


(A. 2) 


are  proportional  to 


M(a,b,z) 


-  U)P  zr 

rJ0  TBT; 


(A. 3) 


where 


(a)r  «  a(a+l)...(a+r-l) 


Since 

M*>  •  /  U-t2)*"-3)/2  dt 

q  -1 

the  results  of  Section  1  show  that 


bq(0)  ■  , 


t 
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Thus 


>,w  ■  «*•!•'> 


(A.5) 


Me  have  already  seen  a  reflection  of  some  of  the  many  recurrence  relations 
satisfied  by  M(atb,z)  but  not  the  useful  differential  property. 


An 

M(a,b,z)  ■  M(a+n  •  b+n  ,  z) 


(A.6) 


dzr 


or  the  Kunner  transformations 


M(a.b.z)  ■  ez  M(b-a  ,  b  ,  -z) 


M(l+a-b  ,  2-b  ,  z)  •  ez  M(l-a  ,  2-b  ,  -z) 


(A.7) 
(A. 8) 


The  Kunner  function  M  Is  related  to  the  Bessel  function  Iv  which  arose  In 
Section  2  by  the  formula 


M(a.b.z)  •  ez/2  r(b-a-SJ(z/4)*'b+l1 
-  (2b-2a-l)  (b-2a)_(-l)r 

£  ■  "J I  ~'f  L.T '  "  1  ~  ~  Ii 


4o  - rUi>Tr  •' -  ,b-w*r  (l/2) 


(A. 9) 


For  |z|  large  , 


lira -a  R-l  (a)Jl+a-b)_  __ 

•  TThHt  (4o  r!  -  1  *  0(111  » 


z  a-b  S-l  (b-aJjU-a),  ,s  S 
- i? - 1  «’*0<|I|5» 


(*.») 


uhpru  th*  upper  sign  i»  t»k«n  If  -u/2<arqzOu/2  ,  th*  loupr  othunriu. 
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From  (A. 5)  and  (A.6),  It  follows  that 


B^(«c)  »  b'(K)  /  bq(«) 


£  log  H(»*  .  q/2  ,  k) 


a(3/2,^21/2,4 

Mll/2  ,  q/2  .  k) 


(A. 11) 


From  (A. 11}  and  (A. 7),  we  see  that 


V-o 


mtmmp 


IA.12) 


and  (A. 11)  and  (A. 8)  give  yet  another  form  for  Bq(-ic)  .  (A. 5)  and  (A.9)  show 


-,<«>  -  ^  eK^z  (}) 


-<¥> 


-  (q*2)  ((q-2)/2)  (-l)r 

2  - r-.-r.-prr-- -  I 


r:  iq 


q/2+r 


(k/2)  . 


(A. 13) 


(A.5)  and  (A. 1^1  show  that,  as  «■*•+»  , 


bq(ic)  »  v h  r((q-l)/2)  e*  ic‘(q'1,/2 

S-l  l(q-l)/2)s(l/2)s  #  . 

{  £  - rr* - -S  k  S  ♦  0(ic  5)} 

s»0 


(A. 14) 


while  If  * 


bqfc)  ■  J*  r(lq-l)/2)(-Kr)“% 

R-l  (H)r(-(q-l)/2)r  , 

{  I  — - p - 11  (-«) 

r«0  r* 


+  0(<-R)  } 


IA.15) 


The  expansions  (A.13),  (A. 14),  and  (A. 15)  could  be  used  directly  to  obtain 
the  related  expansions  for  B^(ic)  .  However  the  procedure  used  In  Section  3  Is 
now  seen  to  be  much  simpler  e.g.  It  does  not  require  this  Appendix. 
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